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n ! Abstract 

•^ . Previous applications of the enclosure method with a finite set of observation 

(-H I data to a mathematical model of electrical impedance tomography are based on 

the assumption that the conductivity of the background body is homogeneous and 
known. This paper considers the case when the conductivity is homogeneous and 
unknown. It is shown that, in two dimensions if the domain occupied by the back- 
ground body is enclosed by an ellipse, then it is still possible to extract some infor- 
^ ■ mation about the location of unknown cavities or inclusions embedded in the body 

^D i without knowing the background conductivity provided the Fourier series expan- 

sion of the voltage on the boundary does not contain high frequency parts (band 

f^ ■ limited) and satisfies a non vanishing condition of a quantity involving the Fourier 

coefficients. 
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The aim of this paper is to reconsider previous applications [5, 6] of the enclosure method 
with a finite set of observation data to inverse boundary value problems related to a con- 
tinuum model of electrical impedance tomography [1, 2]. The point is: those applications 
are based on the assumption that the conductivity of the background body is homoge- 
neous and known. However, from a mathematical point of view, the problem whether or 
not one can still extract some information about unknown discontinuity from the finite 
set of observation data without knowing the exact value of the conductivity is quite in- 
teresting. Proofs of some previous known uniqueness results that employ a finite set of 
observation data, for example, [4] for cracks and [3, 9] for inclusions are based on the 
assumption that the conductivity of the background body is known. This is because they 
start with applying the uniqueness of the Cauchy problem for elliptic equations. 

Besides needless to say, we cannot know the exact value of the conductivity of the 
background body. The inaccurate value causes an error on the observation data and 
therefore on the indicator function in the enclosure method. 



In order to describe the problem more precisely let us start with recalling a typical 
application of the enclosure method with a single set of observation data. 
Let f2 be a bounded domain of R^ with Lipschitz boundary. Let D be an open subset 
with Lipschitz boundary of Vt such that D d Vt and fi \ D is connected. Consider a non 
constant solution of the elliptic problem: 

Am = in fi \ D, 

— = ondD. 
ov 

Here v = (z/i, V2) denotes the unit outward normal vector field on d{VL \ D). The Z^ is a 
mathematical model of the union of cavities inside the body. 

In [5] we considered the problem of extracting information about the location and 
shape of D in two dimensions from the observation data that is a single set of Cauchy 
data of u on dVt. Assuming that D is given by the inside of a polygon with an additional 
condition on the diameter, we established an extraction formula of the convex hull of D 
from the data. The method uses a special exponential solution of the Laplace equation. 
The solution takes the form e-'^*e^^-{'^+«'^ ) where r(> 0) and t are parameters; both 
uo and 00-^ are unit vectors and satisfy u ■ u;-*" = 0. The solution divides the space into 
two half planes which have a line {x\x ■ u = t} a,s the common boundary. In one part 
{x\x ■ uj > t} the solution is growing as r — > 00 and in another part {x\x ■ uj < t} 
decaying. Using this solution, we define the so-called indicator function I^^^± (r, t) of the 
independent variable r with parameter t: 

Jan [ ou ov J 

The enclosure method gives us information about the position of half plane x ■ uj > t 
relative to D by checking the asymptotic behaviour of the indicator function as r — > 00. 
For the description of the behaviour we recall the support function hoiuj) = sup^^^, a; • u. 
Moreover we say that u is regular if the set {x\x ■ u = hniuj)} fl dD consists of only one 
point. 

What we established in [5] is: for regular u there exist positive constants A and 
/i(> 1/2) such that, as r — > 00 

\I.,.4r,0)\r^-e^'-^-^ (1.2) 

provided 

diamD <dis(L',(9fi). (1.3) 

This fact is the core of the enclosure method. Since we have the trivial identity 

4,c.x(r,t) = e-^%,^x(r,0), 

from (1.2) one could conclude that: if t > h£,{uj), then the indicator function is decaying 
exponentially; if t = hoiuj), then the indicator function is decaying truly algebraically; if 



t < hoiijj), then the indicator function is growing exponentially. Moreover from (1.2), we 
immediately obtain also the one line formula 

,. log|4,<^x(r,0)| 

lim ■ = hoiuj). 

T >00 J- 

However this is the case when the background conductivity is known. 

Consider the case when the background conductivity is given by a positive constant 
7. In this case the indicator function should be replaced with 

4,.x(r, t) = e--* J^^ |-^|-e-(-+-")« + ^^e^^^-^^^^A ds. 

Needless to say we obtain the same result as above if 7 is known. However, if 7 is unknown, 
then the term 

e ""* / 7^e"^-('^+*^^)Mrfs 



Ian du 
becomes unknown and therefore one can use only the term 

e-* / ^^e™-("+'"^)rfs (1.4) 

Jon ov 

if M = / on dVt is given. 

The purpose of this paper is to give a remark on the problem: can one still extract 
information about the location and shape of D from the quantity (1.4) in the case when 
/ is given? 

In this paper we show that, in two dimensions if the domain occupied by the back- 
ground body is enclosed by an ellipse, then it is still possible to extract some information 
about the location of unknown cavities or inclusions embedded in the body without know- 
ing the background conductivity provided the Fourier series expansion of the voltage on 
the boundary does not contain high frequency parts (band limited) and satisfies a non 
vanishing condition of a quantity involving the Fourier coefficients. 

2 Extraction formulae 

Let VL be the domain enclosed by an ellipse. By choosing a suitable system of orthogonal 
coordinates one can write 

n ^ {(.,.,) I (g^ (!)%!} 

where a > 6 > 0. In what follows we always use this coordinates system. 

Given uj = (001,002) G S^ set oj-^ = {0U2, — ^i)- Then x- {ou + ioj-^) = {xi — ix2){ooi +ic^2)- 
Let V = e™-(^+^'^^). 

2.1 Preliminary computation 

In this subsection first given f = u\q^ we study the asymptotic behaviour of the integral 

du , 
'j—-vas. 
an Oh' 



However, integration by parts yields 



'y—-vds = 7 / u—-ds — 7 / u—-ds 
an Oh' Jan ou Jan ov 



(2.1) 



and we have already studied the asymptotic behaviour of the second term as described in 
Introduction (see (1.2)). Therefore it suffices to study that of the first term. Since 



dv f 

. u—-ds = t(uji + iLJ2) / uvivi — iv2)ds, 
an ou Jan 

we compute the integral in the right hand side. 

Write 

1 00 

f{9) = f {a cos 9, bsinO) = -ao + X! (^mcosm^ + Pm sin m9) 



(2.2) 



m=l 



where 



1 /■2'r 



1 /■2'' 



ar. 



— / f (a COS 9, b sin 9) COS mOdO, Pm = — / f(acos9,bsin9)sinm9d9. 

TT Jo TT Jo 



Define 



70 = "0/2, 7m = (ttm - ^/9m)/2, 7-m = lm,rn>l. 

Lemma 2.1. VFe /iawe.' if a = b, then 

{ar(cji + icJ2)} 



/ Mf (z/i — iv2)ds = 2Txa^ y^ 
Jan ±^r 



m=0 



ml 



7m+l) 



(2.3) 



if a > b, then 



an 



UV [Ui — 11^2 



)ds = 2'Kab Y, i""Jm{-iVa^ - feV(cJi + iuj2))C^if) (2.4) 



m=0 



where Co(/) = A^'ji + A+71, /or m = 1, 2, 



C„(/) = (A_7„,_i + A+7„,+i) I 4/ + (A„7„+i + A+7„_i) 



and 

Proof. Set 2; = e*^. Since 

i/(acos^, 6 sin I 



a — b 



1 ^1 , 1 
^ 2 Va 6 



a + 6 



cos^ sin^ 



/ cos 9\ f sin9 



2 V a 



and 



c/s = ab^ 



fcos9\ /'sin9^ 



\{ 



+ 



d9. 



we have 

dz 
{ui — iu2)ds = ab{A_z + A+z'-^)^. 

Note also that 

f(a cos 6*, h sin 6^) = ^ Tm-z™ 

and 

Xi — ix'}, = B_z + -B+z^"*" 

where 

^ 2 

Using those expressions, we can write 



ah 

I ^ J\z\ 

Define 



/ uv (ui — iv2)ds 

JdCl 

^7m / _ {A_z + A+z-^)z'^-^exp{T{B_z + B+z-^){L0i+iuj2)}dz. 



Ii{t)= [ z^exp\T{B^z + B+z ^){lui + iLU2)} dz. 

J\z\=l '- ' 

Consider the case when a > b. Using the generating function of the Bessel functions, 
we have 



exp |r(i?_2; + 5+2; ^){uji + iuj2)j = ^Jn y-iVa'^ - h'^T{uJi + 1002)) M/ _ , ;. I ^" 



la — h 
a + b 



and therefore 

i+i 



Ii{t) = 27ni-iy+\Ji+i (^-t./^^^^ricu, + iu;2)) l-tJ 



a + b 



^ a — b 
li a = b, then 

h[r) = 0, / < -2; Ii[t) = 2ni , / > -1. 

Since 

/ uv{ui -W2)ds = — V7„(A_/„(r) + A+/„_2(^)), 
Jan I ^ 

we obtain the desired conclusion. □ 



2.2 Main result 

We denote by E{Q) the set of all points on the segment that connects the focal points 
(— Va^ — b"^, 0) and (\/a^ — 6^, 0) of il. It is easy to see that the support function of the 
set E{Q) is given by the formula hE(n){uj) = \/d?' — 6^|a;i|. 

We say that a function j[ff) = f {a cos 6, b sin 6) of 6 is band limited if there exists a 
natural number A^ > 1 such that, for all m > A^ + 1 the m-th Fourier coefficients «„ and 
Pm of the function vanish. Then we know that Cm{f) = for all m > A^ + 2. 

Now we state the main result of this paper. 

Theorem 2.1. Let ^ be a positive constant. Assume that (1.3) is satisfied. Let uo be 
regular with respect to D. Let f be band limited and u be the solution of (1.1) with u = f 
on dQ. 
(1) Let a > b. Let uo satisfy o^i 7^ 0. Let f satisfy 

00 
Y^{sgnuo,rm^Cm{f)^Q. (2.5) 

m=l 



lim_ - log I /_ -i^^vds\ = max {hoiu;), hE(n){uj)), (2.6) 



The formula 

1 , , r du 
— -log / 7—1 

r >oo J- Jqq (ji/ 

is valid. 

(2) Let a = b. Let f satisfy: for some N > 1 am = Pm = for all m with m > N + 1 

and a% + /9^ 7^ 0. The formula 

1 , . r du 



i f ou 

lim —log! / '-/-;— vdsl = m.ax(hr)(uj),0), (2.7) 

— *°° r Jdn ou 



is valid. 



• We say that a D is behind the line x-uj = t from the direction uj if the D is contained 
in the half plane x ■ uj < t. One important consequence of the formula (2.6) is: one 
can know whether the unknown cavity D is behind the line x ■ 00 = hE{n){^) from the 
direction uj, however, in that case one cannot know the line x ■ cu = hoioj) itself from the 
formula. This shows the limit to extract the whole convex hull of D without an additional 
assumption. 

• The assumption that / is band limited is just for a simplicity of the computation 
and can be relaxed. It is possible to apply directly the saddle point method to study the 
asymptotic behaviour of the integrals in Lemma 2.1 for a / that is not band limited. 
Moreover we want to point out that in a practical situation, one cannot produce highly 
oscillatory voltages on the boundary. This is due to the limit of numbers of electrodes 
attached on the boundary of the body. 

• A typical example of a band-limited / that satisfies (2.5) for all u with a;i 7^ is 
the / given by 

f{e) =AcosNe + BsmNe 

where A^ > 1 and A^ + 5^ 7^ 0. See Remark 2.1 below for this explanation. In general we 
have to choose two /s corresponding to whether a;i > or cui < 0. 

Proof of Theorem 2.1. When a = b, the (2.7) is an easy consequence of (1.2), (2.1), (2.2) 
and (2.3). The problem is the case when a > b. We employ the compound asymptotic 



expansion (see page 118 of [8] for the notion of the compound asymptotic expansion) of 
the Bessel function due to Hankel(see (9.09) and 9.3 of page 133 in [8]): 



(2.^ 



Jm{z) ~ 1 


\nzj 
















/ f 

< COS z - 

\ 


mil '^\^f 
2 4/ ^ 




— sin 1 


Z 


rmr 
2 


" 4y 


oo 
s=0 


yA2s+i{m) 



as z — > oo in |arg2;| < vr — 5 for each fixed S G ]0, tt[ where AQ{m) = 1 and, for s = 1, 2 
A^{m) = ^(4m2 - l^){Am^ - 3^) ■ ■ ■ {Am^ - (2s - 1)^). 



First we consider the case when ui > 0. From (2.8) in the case when z = —i\/ a} — 6^r(ti;i + 
iiO'2) we obtain 

/ 1 \ 1/2 , / ArrP' - 1 1 \ 

^4.) - (^) ."(-.)•>-"* (l - ^— + 0(^)) (2.9) 

as r — > oo. Since / is band hmited, one can find A^ > 1 such that, for all m > A^ + 1 
the m-th Fourier coefficients a^. and /3m of / vanish. Then Cm{f) = ior m > N + 2 and 
from (2.4) and (2.9) we obtain 

r / 1 \i/2 , 

/ uv{ui-iu2)ds = 2nab[ e"e~*'^/^ 

Jan \2-KzJ 

(2.10) 

C . 1 s 7V+1 . N+l 1 ^ 

X ^ + ^ U^ ^-(^) + ^^ ^ ^'^™(^) + ^(7^) • 
Here we claim that 

N+l 

E cUf) = 0. (2.11) 

m=0 

It suffices to prove the claim in the case when 

f {a cos 6, b sin 6) = aj cos j6 + (3j sin j6 (2-12) 

for each fixed j = 1, 2, ■ ■ ■ , iV. Since Em=o CUf) = C.-i(/) + C,{f) + C,+i(/) and we 
have 

cAf) = 0, 



c,.,/,^..,.h/^ ^-^ W^ 



we get 



i-i / I r\ i+i" 



|-</)^K--(^))iM\/:^) M^^ 



Since 



•*- + -*-U^ =»• 



we see that the claim (2.11) is vahd. Therefore (2.10) becomes 

/ uvivi - W2)ds = iirabz'^ e^^e"*"/^ Y m^CM) + O(-) . (2.13) 

Jdn \2tizJ Vm^i ^ / 

Set uji + iuj2 = e'^ with -n/2 < ^ < n/2. Then z^/^ = ^[a^ - 52)i/4gi(^-V2)/2^ gj^^^g 
gi2 ^ erhB(„)(o.)girVa2-fe2^2^ f^^^^ (^^2), (2.1), (2.2) and (2.13) we obtain the compound 
asymptotic formula: 

( du , 

/ -f--vds 
Jan ov 

r— N+1 A 

From this we know that the quantity 

exp\^-Tmax{hD{uj),hE(Q){uj))]\ / ^T^vds\ 

is truly algebraic decaying as r — > oo. Note that we have used the lower bound of fi: 
fi > 1/2. Therefore we obtain the formula (2.6). Next consider the case when ui < 0. 

Write R^{t; f) = / /f(z^i — iv2)ds- Then we have R^{t; f) = —R^^{T;f*) where 

Jdn 
f*{x) = f{—x). Since the ?72-th Fourier coefficients of /* are given by (—1)™ times those 

of / and the first component of —uj is positive, we can derive the corresponding result in 

the case when uji < from the result in the case when uji > hj replacing Cm{f) in the 

condition (2.5) with -(-1)™C^(/). D 

Remark 2.1. Fix j = 1, 2, ■ ■ ■ , iV and let / be given by (2.12). Then a direct computation 

similar to the proof of the claim (2.11) yields 

oo 

E m'CUf) = (j - l)'C,-i(/) + (j + irC.+iif) 



m=l 



~j{a' - 62)-(^-i)/2|(a + by^, -{a- by-f-}. 



This yields also 



m=l m=l 

= (-iy^j(a' - b'r^^-'^'iia + by J, - (a - by^}. 

These yield: a / whose Fourier coefficients aj and (3j vanish for all j > iV + 1 with some 
A^ > 1, satisfies the condition (2.5) if and only if 

N 

J2isgnLU^yJia'' - b''p^-'^/^{ia + by^, -{a- by-fj} ^ 0. (2.14) 



It is clear that there are many /s satisfying the condition (2.14). 

Remark 2.2. In (1) the case when o^i = is not treated. In this case 002 = ^1- li uj2 = 1, 

then from (2.8) we have 

/ 1 X 1/2 



where z = —i\jG? — b"^ r(uji +iuj2). Then from (1.2), (2.1), (2.2) and (2.4) the problem 
can be reduced to the study of the asymptotic behaviour of the quantity 

(2.15) 
as r — > oo. This seems very complicated for general r. However, if we choose 

^"^ / = 1,2,---, (2.16) 



Va2 - 62 
then (2.15) becomes 

(-)'-- { E (i . >^ c„(/) . .: E (-1)- (i - >^ a„(/)) 

= ^^^^ — E ^'{Cn.(/) - ^c'„(r)}. 

^^ m=l 

Note that we have used the claim (2.11) for / and /*. Therefore if / satisfies the condition 

oo 

E^'{C™(/)-^C'„(r)}7^0 (2.17) 

m=l 

instead of (2.5), then for r given by (2.16), the formula 

1 , , /■ du 



lim —log! / 7— — frfsl = max (/i£)(c<j), 0), 
I — >oo r Jan ov 



is valid. By replacing / with /*, we know also that: if UJ2 = —1, then the same formula 

is valid provided 

00 

E^'{^-(/)+^C'„(r)}^0 (2.18) 

m=l 

instead of (2.17). From the computation in Remark 2.1 one can sum the conditions (2.17) 
and (2.18) up in the single form: 

N 

E {1 + (-l)^(sgnc.2)4 j(«' - b'p^~'^/'{ia + 6)^7, -{a- 6)^7"} ^ 
i=i 

where N > 1 and chosen in such a way that, for all m > N+1 the m-th Fourier coefficients 
of / vanish. 



2.3 Uniqueness 

As a corollary of Theorem 2.1 we obtain a uniqueness theorem. 

Corollary 2.1. Let '-f be a positive constant. Assume that D satisfies (1.3). 

(1) Let Q be a domain enclosed by an ellipse. Let /+ and /_ be band limited and satisfy 

oo 
m=l 

Let u± be the solution of (1.1) with u± = f± on dQ. Then the Neumann data 'ydu^/du 
and 'ydu^/du on dQ uniquely determine the convex hull of D U E{Q). 

(2) Let Q be a domain enclosed by a circle. Let f be band limited and non constant. 
Let u be the solution of (1.1) with u = f on dfl. Then the Neumann data 'jdu/du uniquely 
determines the convex hull of D U {0}. 

We emphasize that 7 is unknown. This makes the situation difficult definitely. Assume 
that we have two unknowns {D,'~f) = (Di,7i), (D2,72) and solutions Ui and U2 both 
satisfying (1.1) and the boundary condition u = f on dQ. The key point of a standard 
and traditional approach is to prove that if 'yidui/du = '~f2du2/di' on dQ, then ui = U2 in 
a neighbourhood of dQ. If 71 = 72, then the conclusion is true because of the uniqueness 
of the Cauchy problem for the Laplace equation. However, if 7 is unknown, i.e., the 
assumption 71 = 72 is dropped, one can not immediately get the conclusion (note that we 
are considering a finite set of observation data not the full Dirichlet-to-Neumann map). 
Our approach skips this point by using an analytical formula that directly connects the 
data with unknown discontinuity. 

The proof of Corollary 2.1 is based on: given D the set of all directions that are 
not regular with respect to -D is a finite set; the formulae (2.6) are valid for f = f± in 
(1); the formula (2.7) is valid for / in (2). Therefore, for example, in (1) we see that 
the Neumann data uniquely determine the values of m.a:x{h£,{uj),hE(n){^)) which is the 
support function of the convex hull of D U E{Q) at the directions u except for a finite set 
of directions. Since the support function ho and hE{n) are continues on the unit circle 
and so is max {hoi ■ ), hE{n){ ■ ))• ^ density argument yields the desired uniqueness. 
Remark 2.3. If dD is smooth, then (2) of Corollary 2.1 does not hold. Let Q be the unit 
open disc centered at the origin of the coordinates system and for < i? < 1 let D{R) 
be the open disc centered at the origin with the radius R. Let < Ri,R2 < 1- Fix an 
integer m > 1. For each j = 1,2 let Uj be the weak solution of the problem (1.1) with 
D = D{Rj) and the Dirichlet data Uj{r,9)\r=i = cos m6 where {r,6) denotes the usual 
polar coordinates centered at the origin. Then we know that 

Mr, 0) = Y^^(^" + R'rr'^ cos mO, U2{r, 6) = j^^l^" + ^^^~'") cos me. 

This yields 

1 + Rl"'du2 ^ 1 + Rf^dui 

K — - — = mcos mU = 7, — t: — on oil. 

l-Rl"" du 1 - Rj"' du 

Since i?i and R2 are arbitrary chosen, this means that one cannot uniquely determine 
D{R) from the single set of the Dirichlet and Neumann data f{6) = cos md and 'ydu/du 
on dQ in the case when 7 = (1 + i?^™')(l — i?^™"). This suggests that the singularity of 
dD is essential for the validity of (2) in Corollary 2.1. 
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3 An application to the inverse conductivity problem 

The idea in the proof of Theorem 2.1 can be apphed to the case when the unknown 
domain D is a model of an inclusion. 

We assume that the conductivity k = k{x) of the body that occupies Q is given 
by k{x) = '-f ii X ^ Q \ D; k{x) = 'y ii x E D. It is assumed that the 7 and 7 are 
positive constants and satisfy 7 7^ 7. The vohage potential u inside the body satisfies the 
equation V ■ kVu = in fi. Given uj = (001,002) G S^ set uj-^ = (cj2, — ^i)- Let r > and 

In [6] we have already proved that if u is not a constant function and D is polygonal 
and satisfies the condition (1.3), then for a given direction 00 that is regular with respect 
to D the formula 

1 . , r I du dv 



lim -log| / i-f--v --f--u]ds\ = hoiuj) 
— *«^ T Jan \ ov ov 



is valid. Note that /c = 7 on dVt and we do not assume that the conductivity 7 of D is 
known. 

Here we propose the same question as that of Introduction. Assume that we do not 
know k in the whole domain. Given a non constant voltage potential f = u\qq^ on dVt is 
it possible to extract some information about the location of D from the corresponding 
current density kdu/dv on dVtl 

The answer is yes in the case when the Vt is enclosed by an ellipse. It starts with 
recalling the equation 

du I" dv _ /" dv 

1^-vds = / -f--uds - (7 - 7) / u--ds. (3.1) 

Recall Key Lemma in [6]: there exist positive constants B and A(> 1/2) such that, as 

T > 00 

«^rfs| ~ 4e"^-(-). (3.2) 

Then from (2.2), (3.1), (3.2) and Lemma 2.1 we see that the completely same statements 
as those in Theorem 2.1, Corollary 2.1 and Remarks 2.1 and 2.2 are valid. 
Remark 3.1. In [7] we employed the difference of the values of the voltage at arbitrary 
fixed two points on the boundary of a general two-dimensional bounded domain Q with 
smooth boundary. More precisely we introduced the operator 

Ak{P,Q)--9^^u{P)-u{Q) 

where P and Q are two arbitrary points on dQ; g satisfies / gds = 0; the u is a solution of 

Jan 
the equation 'V ■k'Vu = in f] and satisfies the Neumann boundary condition kdu/dv = g 

on dD.. 

Given ou = (cui,cu2) G S^ set u^ = {0U2, -^1). Let r > and v = e™-('^+^'^^). What we 

have proved is: ii g = dv/dv on dVt and D is polygonal and satisfies the condition (1.3), 

then for a given direction 00 that is regular with respect to D the formula 

^lirn^ ^ log I {Afe(P, Q) - A,{P, Q)} {g)\ = hoiiv), (3.3) 

11 



is valid. Note that we have used the relationship 

{Ak{P,Q)-A,{P,Q)}{g) = -{Ak/j{P,Q)-MP,Q)}ig). 

If 7 is unknown, then one cannot use the term A^{P,Q){g) in (3.3). However, that 
has the simple form 

A,{P,Q){g) = -{v{P)-v{Q)} 

7 

for g = dv/dv on dVt. Using this form. Proposition 3.1 and Lemma 3.1 in [7], one im- 
mediately gets the following formulae provided D is polygonal and satisfies the condition 
(1.3) and uj is regular with respect to D: 

• if a; is not perpendicular to the line passing through P and Q, then 

^lim^-log|Afc(P,(5)(fi')| = max (/i/)(cj), /i{p,q}(cj; 

• if C(J is perpendicular to the line passing through P and Q, choose, for example, 

then 

lim -log|Afc(P,(5)(5')| =max (/iz)(cj), /i{p,q}(c^)) . 

I — >oo T ^ ' 

4 Conclusion 

We confirmed that: in the case when the background conductivity is homogeneous and 
unknown the enclosure method still works provided: 

• the domain that is occupied by a background body has a simple geometry; 

• the Fourier series expansion of the voltage on the boundary does not contain high fre- 
quency parts (band limited) and satisfies a non vanishing condition of a quantity involving 
the Fourier coefficients. 

However, the method yields a less information about the location and shape of un- 
known cavity or inclusion compared with the case when the conductivity is known. We 
found an explicit obstruction that depends on the geometry of the background body. 
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